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Abstract 

The definition of a double-scaling limit represents an important goal in the development of tensor models. 

We take the first steps towards this goal by extracting and analysing the next-to-leading order contributions, 
^f^ in the 1/N expansion, for the IID tensor models. We show that the radius of convergence of the nlo series 

^H coincides with that of the leading order melonic sector. Meanwhile, the value of the susceptibility exponent, 

^D 7nlo = 3/2, signals a departure from the leading order behaviour. Both pieces of information provide clues 

Cn for a non-trivial double-scaling limit, for which we put forward some precise conjecture. 
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1 Introduction 



Growing evidence is being accumulated for the (Tensorial) Group Field Theory ((T)GFT) formalism [1, 2, 3, 4] 

I— I as a promising overarching framework for a quantum theory of gravity; one that is able to incorporate aspects 

^j^ of several current discrete approaches within a powerful quantum field theory setting. TGFTs are theories of 

' rank-D tensorial fields which generate, in their perturbative expansions, a sum over /^-dimensional cellular 

ij (usually, simplicial) complexes. Their simplest incarnation are tensor models [6, 7], wherein the tensors have 

^ r*| finite index sets of size A^. These were proposed already in the early '90s as an attempt to reproduce, in 3d 

'""' and 4(i, the successes of the matrix model formalism in defining both a controllable sum over topologies and 

I a theory of random discrete geometries with a nice continuum limit (given in 2d by Liouville gravity). Such 

^ tensor models describe discrete geometry is purely combinatorial terms (the natural notion of distance being 

"^j" the graph distance on each cellular complex). Their Feynman amplitudes can thus be understood in terms of 

CO the Regge action for discrete gravity evaluated on equilateral triangulations. Moreover, the perturbative sum 

over Feynman diagrams coincides with the definition of quantum gravity given by the (Euclidean) Dynamical 

triangulations approach [10], after appropriate identification of their respective parameter sets. When one 

jNl enriches the combinatorics of tensor models with the group-theoretic data suggested by Loop Quantum Gravity 

[9], Spin Foam models [16] and simplicial geometry [11], one obtains (Tensorial) Group Field Theories: proper 

field theories, with richer state spaces (with generic states being superpositions of spin networks) and quantum 

•• amplitudes, given by simplicial path integrals and spin foam models. It is these richer field theories, building 

. {^ up on the understanding of quantum geometry obtained in loop quantum gravity, that we believe offer the 

S^ most promising candidates for a complete quantum theory of gravity. Actually, with the appropriate data and 

H constructions [13, 8], TGFTs provide what can be argued to be the best fundamental definition of covariant 

loop quantum gravity dynamics, adapted to a simplicial context. In particular, TGFTs provide loop gravity and 

spin foams, as well as dynamical triangulations, with powerful, analytic field theoretic tools, suited to study of 

non-perturbative physics, the dynamics of many degrees of freedom, and the extraction of effective continuum 

geometry. 

While a main motivation for TGFTs is quantum gravity, this is not their only reason of interest. TGFTs can 
be seen, more generally, as a new class of quantum field theories, posing interesting mathematical challenges, in 
particular from the axiomatic and renormalization theory perspective [4]. At the same time, they define a new 
approach to statistical systems on random lattices, such as spin glasses, dimers, Ising and loop models, and in 
this direction they have already produced interesting results [24, 25, 26, 27]. 

As mentioned, and whatever the perspective, their crucial asset is to provide a new setting in which unsolved 
problems can be tackled with the aid of powerful analytical tools from statistical and quantum field theory. In 
fact, many important results have been obtained in the last few years, confirming such potential. It is not the 
place to review all these results [6, 1, 2, 3]. Beyond model building of 4d gravity models, mainly from the spin 
foam and loop quantum gravity perspective, as well as the associated study of their quantum geometric degrees 



of freedom (see [12, 8, 14] and references therein), work in tensor models includes: i) a detailed understanding 
of the combinatorics and topology of the cellular complexes generated in perturbative expansion, which takes 
advantage of results in combinatorial topology [36] , concerns the absence of extended topological singularities 
[37], as well as the presence of embedded Riemann surfaces [38]; ii) the important identification of a large- 
N expansion for tensor models and topological GFTs [18, 19, 20] (other types of large-N expansion have been 
proposed in [21, 22]); leading then to in) many further results concerning the critical behaviour of various tensor 
models [23, 28] and topological GFTs; and iv) the identification the leading order sector as branched polymers 
[29]. Many more results concern field theory aspects of the formalism, including universality [34, 35], scaling 
behaviour [43], renormalizability [39, 41, 42, 44, 45, 47, 46, 17, 15], Schwinger-Dyson equations [30, 31, 32, 33] 
and quantum and classical symmetries [56, 57], non-perturbative aspects [48, 49]. Finally, ways to extract 
effective continuum physics have been explored, for example [51, 52, 53, 54, 55], culminating in the recent [50]. 

Despite all these recent successes, much remains to be done. In particular, given its crucial role in ensuring 
analytic control over the perturbative expansion of these models, it is important that we improve our under- 
standing of the large- A/" expansion of both tensor models and the more involved TGFTs. The first step is to go 
beyond the leading order in such an expansion, which is by now well understood, with the aim to understand 
the next-to-leading order (and possibly yet more sub-dominant) behaviour. In full TGFTs, subdominant pro- 
cesses may become dominant in certain scenarios, e.g. for particular boundary states, or with different choices 
of weights (that is, within different models). Even in the simplest tensor models, control over sub-dominant 
orders is necessary to be able to define double (and multiple) scaling limits. While for higher-dimensional 
models, one should not expect two parameters to control the full series, the two parameters in the simplest IID 
model should at least allow one to extract a broader subclass of graphs than just the leading order graphs. In 
turn, they should capture better the statistical and topological properties of the sum over complexes and reveal 
new critical behaviour, as has been achieved in matrix models [58, 59]. From the perspective of Dynamical 
Triangulations, the aim is to study analytically the continuum limit for finite Newton's constant. 

In this paper, we study the next-to-leading order in the large-N expansion, focusing on the simplest tensor 
models: the independent identically distributed (iid) tensor model, in any dimension. We consider this the 
necessary first step before tackling more involved tensor models or TGFTs proper. To begin, we provide a brief 
review of such IID model, the combinatorial structures arising in their perturbative expansion, some key tools 
for their analysis, the large-N expansion and the leading (melonic) order. We then move on to present the new 
results of our work. Our first main result is that we identify the graphs contributing to the next-to-leading 
order, starting from their core graphs. We show that they correspond to a precise family of graphs decorated 
by melons, generalizing melonic diagrams with a single 2-dipole insertion. We then show that it is possible 
to use the Schwinger-Dyson equations of the model to obtain a closed expression for the connected 2-point 
function at next-to-leading order, as a function of the same quantity at leading order. From this, one can 
extract the critical behaviour of the free energy for next-to-leading order graphs. We show that the critical 
value of the coupling constant is the same as at leading order, and we identify the new critical (susceptibility) 
exponent. This is our second main result. In the process, we unravel a few more interesting technical properties 
of the combinatorial structures generated by tensor models. Together with most of the technical details and 
proofs of the main results, they can be found in a final appendix. We close with an extended discussion of the 
double scaling limit, explaining the implications of our results for this issue, and putting forward some precise 
conjectures concerning its realization. 

1.1 Random matrices, 2d dynamical triangulations and double scaling limit 

Before moving into our tensor model case, we summarize some key results from matrix models. This should 
serve to clarify some of our motivations as well as providing a template of what one could hope to achieve in 
simple tensor models. 

Consider a matrix model, based on a complex N x N matrix M, defined by the partition function: 

ZN,g = f[dMdM] e-^^5(^'^) , where Sg{M,M) = tr(MMt) + ^tr(M^) + ^tr((Mt)3) ^ ^^^ 

M^ = M and g = XX. A Taylor expansion of the associated free energy per (complex) degree of freedom 
{En^q = {l/N'^) log Zn^q) results in a weighted sum over connected Feynman diagrams Q, dual to triangulations 
of orient able 2-dimensional surfaces: 

EMg = y ^77^ 5'''' '/' N-^'^^ , (2) 

'^ /l^ SYM{g) i' ' W 



where hg is the genus of the surface encoded by the graph Q and SYm(^) is a symmetry factor.^ One can 
organize the graphs according to their topology in a 1/A^-expansion: 

E^,, = J2e,,,N-'\ where E,,, = ^ _i^ 5IV.I/2 , (3) 

and it is clear that only the /i = sector of graphs, that is the spherical triangulations, survives in the large- A^ 
limit. Analyzing the series defined by E^^g^ one finds that it has a finite radius of convergence Qc and leading 
order (non-analytic) behaviour given by: 

E^^g - ao f 1 - — j , where ^ = ~2 ' ^^^ 

where the critical exponent 7 is known as the string susceptibility. As one tunes the coupling constant to its 
critical value, g ^ Qc^ a non-perurbative regime is reached, controlled by those graphs with increasingly large 
numbers of vertices. 

A double scaling limit ensues from the fact that the series in g at all orders in the 1/A^-expansion have the 
same radius of convergence ^c, although different critical behaviours given by: 

/ x(2-7)(l-/i) 

Eh,gr^ah(l--\ . (5) 

Up to issues of stability for the series in 1/N ^ if one tunes: 

N ^ 00 , g ^ gc , with N ll- —\ = n , (6) 

where /^ is a constant, then the series F^^g = N'^Ejy^g ~ ^^ ah /^'^^"^^ includes contributions from all topologies. 
This is the double scaling limit. 

To provide these amplitudes with a gravitational interpretation, one defines bare Newton's {G) and cosmological 

(A) constants through: 

1 /3 

log^'=^ and log^ = -— a^A, (7) 

where a is a length, then one may recast the weights into the following form: 



^|Vd/2 ^2-2^5 = exp 



~a'AJ^2 + ^{2Mo-Af2) 



g-5'G,A,a(A/'o,A/"2) /g\ 



where A/2 = \Vg\ and A/q = \Tg\ are the numbers of triangles and vertices respectively in the triangulation 
represented by Q. The exponent on the right hand side is the Regge action for an equilateral triangulation with 
edge length a and thus, is of the form prescribed by the Euclidean Dynamical Triangulations (EDT) approach 
to 2-dimensional quantum gravity. 

In this gravitational re-phrasing, the large-A" limit corresponds to the limit in which (the bare) Newton's 
constant vanishes: G ^ 0. Tuning the coupling constant g to its critical value, and thus to a regime controlled 
by those graphs with increasingly large numbers of triangles, corresponds to the large- volume limit. However, 
by tuning the edge length to zero simultaneously, one can obtain a continuum limit characterized by surfaces 
with finite macroscopic area. The expectation of the area observable A = (\/3/4) a^ A/2 is: 

(A) ^^^ a^^^log^o,, '^' a'(l-^y\ (9) 



So tuning: 

g^gc, a^O, where a^M-^j = 1/Ai? , (10) 

where A^^ is a renormalized cosmological constant. 

The double scaling limit (6), in this perspective, has the advantage of taking into account all 2d topologies at 
the quantum level, but also of accessing the regime of finite Newton's constant. Indeed, it gives a renormalized 
constant: 

1 „. 1 / --X2-7 



^^ _JlogK=- + 81ogM - J) . (11) 



-•^The genus of an orientable triangulated surface is defined as hg = 1 — (\Vg\ — \Sg\ + \Tg\)/2, where \Vg\, \Sg\ and |J^g| are the 
vertices, edges and faces of Q, respectively. 



2 Tensor model essentials 

We now review the basic definitions and properties of IID tensor models, their 1/A^ expansion, and the mathe- 
matical tools that are used to analyze the combinatorics and topology of their Feynman graphs. 

2.1 IID model 

Consider D -\- 1 complex rank-D tensors: 0^., where i G {0, . . . , I)} is the color of the tensor. Moreover, 
each subscript n^ is actually an abbreviation of the form n^ = (n^^-i, . . . , n^o, ^iD, • • • , ^ii+i), where each 
riij G {!,..., N} for some N. The (D + 1)— colored IID model is defined by the partition function: 



Zjy 



,,j = |[##]e-^(^'^), (12) 



where [d(j)d(l)] is the Gaussian-normalized measure on each of the {D -\- 1)N^ tensor components^ and: 

i=0 rii n i—0 n i=0 

It is dependent on three parameters: {N^ A, A}; the size N and two coupling constants. In the interaction terms, 
J2n denotes the sum over all indices n^j, subject to the condition that riij = riji. Thus, each tensor shares one 
argument pairwise with each of the other D tensors. Let us remark briefiy that this colored IID model, defined 
in terms of D -\- 1 complex tensors with simplicial interaction^, is equivalent to a tensor model for a single 
tensor. This equivalence may be directly constructed via successive integration of all but one the tensors within 
the partition function of the colored simplicial model above and leads to an effective action for the remaining 
tensor that contains an infinite number of [/(A^)^ -invariant interactions, whose respective coupling constants 
are precise monomials of A A (see [35] for details). 

2.2 1/N— expansion 

One recognizes immediately that expressions such as (12) are not naively integrable, so one performs a Tay- 
lor expansion of the integrand with respect to the coupling constants, A and A, to obtain more manageable 
quantities. One evaluates the resulting Gaussian integrals via Wick contraction. The result is summarized in a 
Feynman expansion as: 

g = AA, 

where the Feynman graphs Q label the pattern of contractions, SYm(G) is a symmetry factor, uj{Q) is the degree 
of divergence, while \Vg\ and \Tg\ are the number of vertices and faces in Q, respectively. 

There are two of remarks to be made at this stage: 

— The Feynman graphs Q are closed (D + 1)— colored graphs. This coloring allows to encode topological 

information, in such a way that such graphs are topologically dual to abstract simplicial D-dimensional 
pseudomanifolds. We will give more details on the definition and properties of colored graphs in the 
following. Thus, the partition function is a weighted sum over such objects. 

— The degree is a non-negative graph-dependent integer. Therefore, graphs may be ordered according to their 

degree and since it is bounded from below by zero, it makes sense to consider a 1/N— expansion. 

In fact, one can re-organize the graphs as: 

^N,, = iV^EE^-.^^"^^^V, where Z^,, = ^ -i— . (16) 



^ P g . C^(^) 



SYM(g) 

|Vg|=2p 



^By Gaussian-normalized, we mean: 



•' \ i=0 n / 



, , 1, (13) 

"^The pairing of indices mimics the gluing of {D — l)-simplices across common faces to form a D-simplex. 



Calculating the coefficients Z^^p allows one to extract the critical behaviour of the series: 



(17) 



that is, the behaviour of the partition function at a given order in the l/A/'-expansion. With this in mind, one 
must label and enumerate the graphs at the order of interest, which in turn requires a more detailed examination 
of the {D + l)-colored graphs, to which we now turn. 



2.3 Essentials of (D+l)-colored graphs 

In this section, we present an intuitive description of some basic features of {D + l)-colored graphs. For more 
technically precise definitions, we refer the reader to [6]. 

(D+1)— colored graphs: A (Z^ + l)-colored graph is a graph comprising of {D + l)-valent vertices, such that 
any given vertex is colored either black or white, and each of its D -\-l incident edges is distinctly colored 
from the set {0, ...,!)}. Moreover, the vertices are connected so that black vertices have only white 
neighbours and vice versa. An example is provided in Figure 1. 




Figure 1: A {D -\- l)-colored graph, with I^ = 3. 

Importantly for applications to quantum theories of gravity, these colored graphs are topologically dual 
to D-dimensional abstract simplicial pseudo-manifolds. 

k— bubbles: One identifies the k— bubbles of species {io, . . . , ik-i} as the maximally connected subgraphs 
containing the k distinct colors: {zq, • • • , ik-i} C {0, . . . , D}. In an obvious fashion, /c-bubbles are nested 
within {k + l)-bubbles and so forth. More subtly, the /c-bubbles are dual to the {D — /c)-dimensional 
simplices in the associated simplicial complex, while the nesting relations encode how these simplices are 
glued together. Some /c-bubbles of Figure 1 are identified in Figure 2. 




Figure 2: The 3-bubbles of species {0}, {!}, {2} and {3} (clockwise). 

In particular, note that the faces are the 2-bubbles. We shall often use the notation: {zq, • • • ,^/c-i} 
{0,...,i:>}\{zo,...,^/c-i}. 



It is possible to express the degree of a graph in terms of the degree of its {k -\- l)-bubbles: 



{io...ik;p) 



where {io - - - ik] p) labels distinct (A:+l)-bubbles in Q and ^[^+^1 denotes the total number of (A:+l)-bubbles 
in Q. which in the case of its I^-bubbles, reduces to the following relation: 



^{G) = ^^y^ (p + D- Bl^l) + ^ a;(B(.^^)) . (19) 



k— dipoles: One wishes to catalogue graphs. For colored graphs, the key tool to do so is a class of combinatorial 
moves that have a well-controlled effect on bubble structure. These transformations are known as k— 
dipole moves. A /c-dipole move of species {zq, • • • ,^/c-i} is illustrated in Figure 3. As one can see, there 
are actually two types of dipole moves, dipole creation and dipole annihilation, one being the inverse 
of the other. A /c-dipole annihilation consists, roughly speaking, in the removal of k lines connecting a 
white and a black vertex, together with the vertices themselves, while joining the remaining D^l — k lines. 
There is one condition that must be satisfied by the k edges of colors ii, . . . , z/^ on the right of Figure 3 - 
they should separate two distinct (I) + 1 — /c)-bubbles of species {io, • • • ,i/c-i}- Thus, /c-dipole creation 
(annihilation) increases (resp. decreases) the number of (i^ + 1 — A:)-bubbles of species {zq, • • • ,i/c-i} by 
1. 




Figure 3: The k-dipole moves of species {zq, • • • , i/e-i}- 

Additionally, if the (1^ + 1 — /c) -bubble added (removed) is a (I) — /c) -sphere, then the /c-dipole implements 
a homeomorphism on the associated topological space. 

Then, consider the following example where Q possesses a /c-dipole of some species and the graph resulting 
from the annihilation of this dipole is denoted by Q\dk. Their respective degrees are related by: 

uj{g) = uj{g\dk) + ^^~^^' {k - 1){D - k) . (20) 

Importantly, for both k = 1 and k = D the degree is unchanged. 

1— dipoles and core graph equivalence classes: Given that 1-dipole moves preserve the degree, it is per- 
haps unsurprising that they play a special role in cataloguing Feynman graphs of the IID model. One can 
partition the graphs with a given degree into equivalence classes, where the 1-dipole moves constitute 
the equivalence relation. In other words, two graphs are in the same equivalence class if they are related 
by a sequence of 1-dipole moves (both of creation and annihilation type). Furthermore, it emerges that 
these equivalence classes come equipped with convenient representatives, known as core graphs - those 
members of the class from which no more 1-dipoles can be annihilated. In general, there are several core 
graphs within a particular equivalence class. However, this does not pose a problem. One simply picks 
one such graph for each equivalence class. The rest of the graphs in the class are generated by performing 
arbitrary sequences of 1-dipole moves on this core graph. 

One would also like to label each graph in the equivalence class uniquely, so that the coefficients in (16) can 
be computed. This turns out to be tricky and, unfortunately, the sequences of 1-dipole moves mentioned 
a moment ago are not the best tool to achieve this goal. The reason is that often there are several distinct 
sequences that transform a representative core graph to the same graph in the equivalence class. The next 
section details the solution to this problem in the leading order and next-to-leading order sectors. 

Jackets: Jackets are the name given to a certain class of 2-dimensional surfaces embedded within the D- 
dimensional topological manifold. They are encoded via a (I) + l)-cycle a of the set {0, . . . ,D}. The 
surface is constructed from the cycle as follows. Consider a graph Q and a planar projection of the 
neighbourhood of each black vertex such that the incident colored edges are ordered clockwise around the 
vertex according to the cycles a. For white vertices, the cycle determines the anti-clockwise ordering. 
Since a (I) + l)-colored graph is not generally planar, away from the vertices the edges cross and the 



graph is embedded in a surface of non-zero genus. This surface is provided by Q itself, comprising of the 
totahty of its vertices and edges, along with all the faces whose two colors are adjacent in the cycle a. 

There are D\/2 distinct jackets J" in a graph. ^ As a result, the degree may be re-expressed as a sum over 
the genera of the jackets:^ 

uj{g) = Y.hj. (21) 

J 

Gravitational interpretation: As in the 2d case outlined in the introduction, the amplitudes have a gravita- 
tional interpretation in terms of the Regge action evaluated on equilateral triangulations^ , and the sum 
over graphs can thus be put in correspondence with the definition of quantum gravity suggested by the 
Euclidean Dynamical Triangulations approach. Defining: 



one may recast the weights as: 



-VdAA/'d + -—^ (27rVD_2A/'D-2 -^ ^Vd_2 A/}) arccos — 

IottG \ 2 D 



(23) 

where JVd = \Vg\ and Md-i = \^g\ are the numbers of D- and {D — 2)-simplices respectively in the 
triangulation represented by ^, while Vk = {a^ /k\)^{k + l)/2^ is the volume of an equilateral /c-simplex 
with edge-length a. Once again, this is the action prescribed by the EDT approach. The large-A^ limit 
corresponds to the vanishing of (the bare) Newton's constant: G ^ 0, while tuning the coupling constant 
to its critical value lead to a regime whose behaviour is controlled by D-complexes with increasingly large 
numbers of I)-simplices. As in the 2d case, a double scaling limit would then allow not only to include 
a more general class of triangulations in the sum (although in this case it may not allow to go beyond 
spherical topology), but also to probe the regime corresponding to finite Newton's constant. 

2.4 Observables 

Rather than deal with the partition function directly, two other observables are studied in this paper. First, 
the free energy is defined as: 

EN,g = ^\0gZN^g, (24) 

and its contributions come from connected closed {D + l)-colored graphs. Meanwhile, the connected 2-point 
function is defined as: 

{<i>in C)c = -^ l[dcf>d^] cf>i^ 4^^ e-^(^>^) , (25) 

^N,g J 

and its contributions come from connected {D + l)-colored graphs with two external edges of color i. Since all 
connected closed {D + l)-colored graphs are also 1-particle irreducible (iPi), cutting a single edge of color i 
within any connected closed graph gives a connected 2-point graph. Moreover, for a closed graph Q scaling like 
TV - {D-iy.^y^) ^ its associated 2-point graph, ^, scales like N~ ^^-^y-^^ K In other words, all graphs get rescaled 
by the same factor N~^ . Thus, a 2-point graph contributes at a certain order to the 2-point function if and 
only if its associated closed graph contributes at that order to the free energy. 

Due to index conservation, the connected 2-point function may be factorised as: 

{(^ln4>^fh)c = GN,gSmrn , (26) 



^Although there are D\ distinct cycles for a set with D -\- 1 elements, reversing the cycle does not does not generate a different 
surface. 

^For an orientable surface: hj = 1 — {\Vj\ — \Sj\ + \J^j\)/2^ where \Vj\^ \^j\-> \^j\ ^-^e respectively the vertices, edges and 
faces of the jacket J . 

^In proper TGFTs, on the other hand, thanks to their richer set of data, the correspondence can be improved to give generic 
simplicial path integrals for discrete (1st order) gravity actions, with generic assignment of geometric variables (areas of triangles, 
holonomies of discrete gravity connections, etc), in turn dual to spin foam models [1, 2, 13, 12, 8]. 



where the factor Gn.q is independent of the color of the external edges. Indeed, there exists a Schwinger-Dyson 
equation relating these observables:^ 

GN^g = l^9^EN^g. (27) 

Thus, the behaviour of the free energy is directly and easily related to that of the connected 2-point function. 
This is a very useful property for the analysis to be detailed below, since connected 2-point graphs are more 
easily catalogued than closed graphs. 

3 Graphs 

We now present the analysis of the combinatorial structure of the graphs appearing in the perturbative expansion 
15, in the 1/A^ expansion, at both leading and next-to-leading order. We focus on the main steps of the analysis 
and on the results, leaving the detailed proofs to the appendix. 

3.1 Leading order 

We shall just state the results obtained in [23]. One finds that the leading order in the l/A/'-expansion is 
specified by: ijo{Q) = 0. Thus, the pertinent coefficients are the: 

G'lo,p '= Go,p = y. 7^' (^^) 

where Q is the closed graph obtained from the 2-point graph Q by joining its two external lines of color 0. 

For the leading order closed graphs, there is a single equivalence class with a unique core graph called the 
supermelon. It is illustrated in Figure 4. 




Figure 4: The supermelon core graph, denoted Q^ 



SUPERMELON- 



Then, one turns to the 2-point graphs. The graphs occurring at leading order are known as rooted melonic 
graphs. The fundamental building blocks of any rooted melonic graph are the elementary melons, illustrated 
in Figure 5. 




Figure 5: An elementary melon of species zq- 

Such a melon consists of two vertices sharing D edges. Both vertices have one external edge. Obviously, both 
external edges possess the same color, say i^. Thus, one refers to such an object as an elementary melon of 
color io. Moreover, an elementary melon has two distinguished features: i) an external edge of color io incident 
to the white vertex, which is known as the inactive edge; ii) D -\- 1 edges incident at the black vertex, which 
are known as active edges. 

The set of all rooted melonic graphs, denoted by M, is the union of the subsets Mp containing rooted melonic 
graphs with 2p vertices: 

M=[JMp. (29) 

p>i 



^The appropriate Schwinger-Dyson equation is: 



One may define the elements of Mp as follows: 

p = 1: There are only 1^ + 1 rooted melonic graphs in Mi, the elementary melons illustrated in Figure 5 for 
different choices of io ^ {0, . . . , D}. 






Figure 6: An elementary melon of color 2 inserted along the active edge of color 2 (for D = 3). The active 
edges are drawn using full lines. 

p = 2: One obtains the graphs in M2 from the graphs in Mi by replacing an active edge of a given color by an 
elementary melon of the same color, as shown in Figure 6. (In fact, this is I)-dipole creation.) 

p = k: One obtains the graphs in Mp from those in Mp_i by replacing some active edge by an elementary 
melon. 

For a graph occurring in Mp, the initial combinatorial factor coming from the Taylor expansion is 1/p!, while 
the graph is obtained from exactly p\ Wick contractions. Thus, the final combinatorial factor is SYm(^) = 1. 
As a result, the problem of calculating the coefficients Glo,p has been reduced to the enumeration of distinct 
patterns of melonic insertions. 

Remark: A generic leading order 2-point graph is obtained from a leading order 2-point graph at p = 1 
by performing some sequence of 1-dipole moves, just as a generic leading order closed graph is obtained from 
SsuPERMELON by somc sequence of 1-dipole moves. However, in our definition of rooted melonic graphs, we have 
only talked about inserting melons within melons. This stems from the result, proven in [23], that: The set of 
rooted melonic graphs M is closed under 1-dipole creation and annihilation. 

3.2 Next-to-leading order: statement of the results 

In this section, we present the results of our analysis of the next-to-leading order sector. One can anticipate its 
components: i) the identification of NLO core graphs; ii) an iterative procedure to generate all graphs at that 
order, starting from the core graphs. For the technical aspects, we refer the reader to Appendix A where the 
precise statements are laid out and proven. 

Our first main result concerns the core graphs: 

Proposition 3.1. Consider the IID tensor model with D > 3. The graphs S2 dipole (seen in Figure 7) are the 
NLO core graphs. 




Figure 7: The core graphs at NLO, denoted collectively by ^2 dipole- 

In the NLO sector, there are { 2 ) ^^^^ graphs, all of the form given in Figure 7. Note that they may be 
obtained from the supermelon graph by creating a single 2-dipole, illustrated in Figure 8, of which there are 
{^2^) distinct species. 

It emerges that the creation of the first 1-dipole in the NLO core graphs is equivalent to the insertion of an 
elementary melon. However, the creation of a second 1-dipole has two possible effects: i) it may again be 
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Figure 8: A 2-dipole of species {ioii}- 



equivalent to the insertion of an elementary melon or ii) it may produce a graph of the form illustrated in 
Figure 9. 





Figure 9: The creation of a second 1-dipole, of species {^o}, which is not equivalent to the insertion of an 
elementary melon. The edges contributing to the 1-dipole are tagged. 

One may iterate this procedure to arrive at graphs of the form drawn in Figure 10. Note that the graphs 
produced by £ — 1 iterations of this procedure have two faces of species {ioh}^ each with 2£ edges. We shall 
denote such graphs by ^2 dipole,^ , in which case, ^2 dipole,i = ^2 dipole- 




Figure 10: The graphs ^2 dipole,^- 

Definition 3.2. One denotes by Si^n^ the set of closed graphs derived from ^2 dipole,^ by inserting arbitrary 
combinations of n elementary melons. 



Proposition 3.3. The set of graphs: 



S=[j Si, 



(30) 



£>1 
n>0 



is closed under 1-dipole creation and annihilation. 



As at leading order, it is easier to accurately count distinct NLO connected 2-point graphs rather than closed 
graphs. Given that the next-to-leading order sector in the l/A/'-expansion is specified by: uj{Q) = cc;(^2 dipole) = 



^^^{D - 2), one finds that: 



G 



NLO,p 



G (D-l)\ 



{D-2),p 



E 

~ 2 

\Vg\=2p 



1 



,n TV SYM(g) 

u;(a)=(^^(D-2) ^^^ 

y : 



(31) 



There are a number elementary building blocks that are used to define 2-point graphs at this order. Supple- 
menting the elementary melons, there are the 2-point insertions obtained by cutting a edge of ^2 dipole,^ (with 
^ > 1). Depending on the edge cut, they take the forms illustrated in Figure 11 and we shall refer to them all as 
elementary 2— dipoles. Note that, once again, the solid edges are active, while there is still just one inactive 
edge, marked by a dashed line. 
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Figure 11: The elementary 2-dipole insertions 7i,o,o- 
Proposition 3.4. For the IID model, the set of NLO connected 2-point graphs is: 

£>1 
m>0 
n>0 



The subsets T^ ^nn are defined as follows: 



771 = 0, n = 0: The graphs in T^,o,o are the 2-point graphs obtained by cutting the edges of ^2 dipole,^, that 
is, they are the elementary 2-dipoles illustrated in Figure 11. 

771, n = 0: The graphs in T^,rn,o are obtained by replacing an interior active edge of an elementary melon 
with a graph from T^^^_i^o- Thus, they have the generic form drawn in Figure 12. 





Figure 12: The recursive definition of the elements in T^,rn,05 where 7£,m,o ^ ^^,m,o and 7^,rn-i,o ^ ^^,m-i,o- 

% m n: The graphs in T£^rn,n are obtained from those in T£^rn,n-i by replacing an active edge with an elementary 
melon. 



Thus, the subset of graphs with 2p vertices is: 



[j T£^rri,p-2£-m 

l<£<lp/2\ 
0<m<p-2£ 



(32) 



As before, the initial combinatorial factor coming from the Taylor expansion for a graph in Si^rn,n is 1/(2^ + 
771 + 7i)!, while the graph is obtained from exactly {2i -\- m -\- n)\ Wick contractions. Thus, the final combinatorial 
factor is SYm(^) = 1. This concludes the identification of the graphs contributing to the next-to-leading order. 

Remark: We stress here the important point that the core graphs S2 dipole and thus the whole NLO sector 
correspond to /^-dimensional cellular complexes of spherical topology, just like the leading-order graphs. As a 
result, the 2-point graphs represent the D-dimensional ball. 



4 Critical behavior 



One now turns to an analysis of Ejy^g, which as one may recall is the free energy at given values of the expansion 
parameter N and coupling constant g. This may be expanded in both parameters: 



'N,g 



^^E^^pN-'^gP where 



OJ p 



Euj,p = Yl 



1 



G: 



\Vg\=2p 



SYM(g) 



(33) 
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The quantity E^^p counts the number of closed {D + l)-colored graphs with a given degree and a given number 
of vertices (weighted by the relevant symmetry factors). In turn, the large-p behaviour of E^^p provides the 
radius of convergence of the series along with the critical exponent:^ 

E^,g = ^E^^pgP^A^(l-^^] as g ^ gc,o. . (34) 

p \ y / 

where A^ is a constant of proportionality, while gc^u and j^ are the radius of convergence and the susceptibility 
exponent, respectively. In the coming section, both leading order and next-to-leading order sectors are analysed. 
The leading order melonic sector not only provides an invaluable introduction to the techniques used, but also 
some necessary results. So, it is worth reviewing explicitly, albeit briefly. 

As already mentioned, the strategy involves examining the behaviour of the connected 2-point function at the 
relevant order. The corresponding behaviour for the free energy can then be found by integrating equation (27). 
To succeed, one needs also the 1-particle irreducible (iPi) 2-point function, given by: 

( ^lui ^m,i )lPi = ^mn ^N,g , (35) 

where ^N,g is a constant, depending on N and g. There is a convenient identity relating this to the connected 
2-point function: 

GN,g = {l-^N,gy^ • (36) 

This leads immediately to the following relations:^ 

^NLO,g — [^ ~ ^LO,g) ^NLO,g [^ ~ ^LO,g) — ^LO,g ^NLO,g ^LO,g 

One simply needs to find some more equations to close the system and solve it for the desired function. The 
equation to be used depend on the detailed combinatorial structure of the graphs at each order. 

4.1 Leading order sector 

The second equation for the leading order 2-point functions descends directly from the melonic structure of the 
contributing graphs. A rooted melonic graph has the generic form given in Figure 13. 




Figure 13: A generic LO 2-point graph. 

Thus, any leading order iPl 2-point graph has the form given in Figure 14, where the shaded circles indicate 
the insertion of an arbitrary connected rooted melonic graph. 

The diagram illustrates the following mathematical relation: 

5^Lo,s=5Gfo,g. (39) 

Using (38), this leads to a closed equation for Clo,^- 

G,o,g = l+ffGfo+^ (40) 

which in turn can be solved via a series expansion for the coefficient: 



^For clarity, we have assumed that there are no logarithmic factors contributing to the leading divergence as ^f — )- gc,uj- In 
practice, one should demonstrate this explicitly. 
^We have used the shorthand: 



X„^o,g = X(£^^i^_^^^g 



(37) 



where X may be replaced by the suitable observable e.g. {E, G, E}. 
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Figure 14: The generic graph contributing to the LO iPl 2-point function. 

where Cp are the {D + 1)-Catalan numbers. Applying Stirhng's formula^^ to the series coefficients, one 

can examine their large order behavior: 



(43) 



Glo,p ~ /?Lo (ffcLo)"^ P~^ where < 

(D + 1)^+1 
Such a series has a radius of convergence gc,Lo and the behaviour of the series in the vicinity of (/c,lo is given by: 



e lD + 1 

Plo 



D 

9c,LO 



i/c,LO 

This implies the following critical behaviour for the free energy: 



GLo,,-^(l-r^) . (44) 



/ Q \^"^^° 1 

E^o,g -- Ao 1 ^— where 7lo = 7^ • (45) 

\ 9c,LO J ^ 

In this context, 7lo is known as the entropy exponent or susceptibility. 

4.2 Next— to— leading order sector 

The analysis at NLO proceeds similarly. In effect, it requires one to obtain the cardinality of the set T defined 
in Section 3.2. Consider a iPl 2-point graph occurring at NLO. It receives contributions from 2-point graphs of 
the form drawn in Figure 15. The first corresponds to the insertion an elementary 2-dipole into an elementary 
melon followed by melonic insertions thereafter (the graphs in Ti^rn.n with m > 1). The others correspond to 
melonic insertions into the elementary 2-dipole graphs (the graphs in T^,o,n)- 

In terms of generating functions, this relation translates into the following equation: 

S..O,. =9D [G.o,.]^-^ C..O,. + ^^^ E^'' [G.oJ''''^''-' , (46) 

where the D-dependent factors arise from the choice of species. After rearrangement, one finds: 

2 2^9 [^hO,g\ 2 9 [^hO,g\ 2^ l^LO,^ - -Lj 

where the second equality uses equation (40). The summation over i can be performed explicitly,^^ since 
|Glo,^ — 1| < 1 in the range < ^ < ^c,lo and one gets: 

DiD±l)_2 if^ i2D+2^^ ^x-1 

^NLo,^ - 7 ^^-^x • (49; 



[l-gD[G,o,gf^') 



^Stirling's formula states that: 



lim , \ ,„ = 1 (42) 

n^o^A/2^(^) 

-'^-'^ Inverting equation (40) for ^f as a function of Glo, one finds: 

3(Glo) = (G',o-1)/[G'.o]''+' . (48) 

One can see that g = corresponds to Glo = 1- Moreover, one knows that the right hand side is a monotonically increasing 
function of Glo up to some critical value Glo, critical 5 where one encounters a stationary point. By examining its derivatives, one 
finds that g (Glo) has a maximum at Glq, critical = {D -\- 1)/D, leading to g (Glq, critical) = D^ /{D + 1)^+-*^ = 5'c,lo, as expected. 
Since | Glq, critical ~ 1| = 1/^ < I5 we have what we need to resum the series in (47). 
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E 




+ E E 

£>1 zi : iij^i 





Figure 15: The NLO consistency equation. 



Moveover, upon differentiating (40), one finds: 

d 



dg 



[G.or^' 



so that: 



^NLO,g — 9 



D f 1 \ d 



2 VGLo,g-2y dg 



[[G.o,,f^') 



(50) 



(51) 



Knowing the critical behaviour of Clo,^ ahows one to determine the behaviour of the NLO series from (51): 



^NLO,5^ ~ I 1 



9c,LO 



since Glo q ^ constant + ( 1 



This imphes the fohowing critical behaviour for the free energy: 

^ "Tnlo 



-t^NLO,g ^ ^NLO I 1 
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9c,LO 



where 7nlo = 



(52) 



(53) 



We find then that the critical point of the NLO is the same as that of the leading order, while the critical 
exponent differs. This is exactly the property indicating the potential for a double scaling limit. 

5 Discussion: a double scaling limit? 

Before concluding, we would like to discuss briefly the implications of our results for the existence and nature 
of the double scaling limit. With the above analysis, we have seen that: 



EN,g - Ao N^ ( 1 



• Ant^o N- 



2-D 



9c 



(54) 



Obviously, these two pieces of information already allow one to conjecture a possible form of double scaling 
limit, whose actual realization relies on the following two properties: 
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Al: There is subset (perhaps with infinity cardinahty) of orders, whose elements are labehed by m G No, 
behaving as: 

E^,g ~ Ara iV^^^-i^) (l - - y (55) 

A2: Ah other orders are stih washed away in the double scaling limit. 
If one assumes these two properties, then the double scaling limit of (A^^^^~^^/^£^Ar,c/) ensues from: 

N ^ oo, g ^ Qc, such that TV I 1 - — j = n . (56) 

We do not have a proof for the above assumptions. The proposed form of double scaling remains, therefore, a 
conjecture. We can give, however, some supporting arguments for it. 

As regards (Al), there are certainly core graphs weighted by N^^'^~-^\ namely, those that reduce to the 
supermelon graph through the annihilation of m successive 2-dipoles.^^ 

Having said that, one would need to show that the sectors, which we call 2— dipole sectors, generated from 
these core graphs each give rise to a resummable series with the radius of convergence and critical exponent 
conjectured in (55). 

Moreover, for a given tti, there may be other sectors of graphs weighted by A^^(2~^). If such sectors exist, one 
would need to show the realization of one of the following two scenarios, either i) they have the same behaviour 
as the 2-dipole sectors or ii) they do not have the same behaviour, but the radius of convergence is larger then 
Qc^ so that they are washed away in the limit. 

To prove (A2), one would need to show that the 2-dipole sectors dominate over all other sectors. This appears 
a rather daunting task. One might attempt to tackle it by showing that for sectors generated from core graphs 
with, say, 2p vertices, then the 2-dipole sector generated from core graphs with 2p vertices dominate. 




Figure 16: A subdominant core graph with four vertices. 

For example, we can show that this is already for those sectors generated from core graphs with four vertices. 
As illustrated in Figure 16, all such core graphs are obtained from the supermelon core graph after the creation 
of a /c-dipole of some species with k >?> and denoted here by S/c dipoie- The techniques that we developed for 
the NLO sector work on these sectors perfectly. ^^ One finds that while their weight in the l/A'-expansion is 
-j\jD+{k-i){k-D) ^ they generate a resummable series with behaviour (1 — g/gcY^'^ . Thus, their contribution to 
the series [N^^^-'^^^'^EN.g) is: 

^3(D-2)/2+(/c-i)(/c-D) (^1 - A y ^ (57) 



in the double scaling limit defined by (56). The aim is to extend such arguments to higher orders in the 
expansion. 

There are two points of further interest pertaining to the double scaling limit conjectured above. 

The first is a comparison to the double scaling limit of the single matrix model. In that context, the double 
scaling limit captures contributions from all topologies. For D = 3, rather than contain contributions from 
all topologies, the limit above captures the complete spherical sector of the IID tensor model, since graphs 
representing the 3-sphere may be obtained from ^supermelon by some sequence of 1-dipole and 2-dipole moves . 
For D > 3, the limit captures less and less of the spherical sector, since performing /c-dipole moves (with k > 2) 



■"^^A subtlety of 2-dipole annihilation is that one may need to insert 1-dipoles briefly in order to effect the annihilation. 

-"^^In fact, for the sectors generated from ^fc-dipole with /c > 3, the equivalent of the closed graphs ^2 dipole, € with -^ > 1 do not exist. 
Thus, the set of closed graphs for these sectors is simply S = U^>o»S'i,n where Si^n = {Cfe-dipole with arbitrary melonic insertions} . 
This, however, is already enough to determine the same critical behaviour. 
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on ^suPERMELON ^lay also result in a graph encoding the I^-sphere. This means that, the higher the dimension, 
the harder it is for the simplest tensor models to capture the topology of the continuum spacetime they aim to 
describe, not to mention its geometry. Richer models involving more coupling constants, thus multiple-scaling 
limits, or more structured amplitudes, depending on a richer set of data, are then called for. 

The second utilizes the gravitational interpretation provided by Euclidean Dynamical Tri angulations. As before, 
it allows one to enter a regime of finite Newton's constant, by defining a renormalized constant: 

1 "'"S" ' . « .,„,A_£V"""'. ,58) 



Gr Vd-2 G Vd-2 V 9^ 

However, for D > 2, Newton's constant is dimensionful. As a result, this Gr is finite in the large-volume limit 
rather than the continuum limit, where one has also a ^ and which will require a more subtle analysis. 

6 Conclusions 

We have studied the next-to-leading order in the large-N expansion of IID tensor models, in any dimension. 
We have identified the graphs contributing to it, corresponding to families S^^n of graphs. We then studied 
the critical behaviour of the free energy for such graphs. The result is that the critical value of the coupling 
constant is the same as at leading order, while the critical (susceptibility) exponent is different. These results 
support the possibility of a double scaling limit capturing more properties of the sum over triangulations for 
spherical topologies, and suggest the form that this double scaling may take, as we have discussed in some 
detail. The analysis we performed for such simple tensor models can also form the basis of a similar analysis 
for more involved tensor models as well as for proper TGFTs. 

A Next-to-leading order: technical details 

In this appendix, we collect the needed technical definitions, report the details of our results and present the 
corresponding proofs. 

A.l Some features of (rooted) melonic graphs 
A. 1.1 Irreducibility 

To begin, we generalize the idea of n-particle reducibility to include some color information. 

Definition A.l. A graph is said to be (ao, . . . , ai))-particle irreducible, if it remains connected, having cut aj 
edges of color j, for all j. 

We shall denote by ej^...j^ the vector with components ji, . . . , j/e equalling and with the rest equalling 1. Here 
are some easily verified facts about {D + l)-colored graphs: 

Ql: All connected closed {D + l)-colored graphs are e^-particle irreducible, for all j G {0, . . . ^D}. In other 
words, one may cut along any D distinctly colored edges and the graph remains connected. 

Q2: All iPl 2-point {D + l)-colored graphs with external edges of color {z}, other than the elementary melon, 
are e^-particle irreducible, for all j G {i}. The elementary melon is the exception, since it does not have 
an internal edge of color {i}. Rather the elementary melon is of species {i} is e^j-particle irreducible, for 
ah JG{?}. 

Now we specialize to melonic graphs. 

A. 1.2 Melonic vertex pairs 

Any closed {D + l)-colored graph has equal numbers of black and white vertices. Thus, in principle, there are 
many possible ways to partition these vertices into black-white pairs. However, certain graph properties serve 
to distinguish particular pairings. 
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Definition A. 2. Rooted melonic graphs have a melonic vertex pairing defined at: 

p = 1: An elementary melon has just two vertices and these form a melonic vertex pair. 

p = k: Such a graph is constructed via the iterative insertion of elementary melons. Each elementary melon has 
a pair of vertices. Thus, the vertices of the graph are paired according to the elementary melon, within 
which they were inserted. 

As a result, for a given vertex in a rooted melonic graph, one may identify its paired vertex by deleting all 
sub- nested melons. Moreover, for a closed melonic graph Q, one may construct a melonic vertex pairing by 
inserting a fictitious cut along one edge and utilizing the resultant rooted graph. From this, one can derive a 
number of properties: 

PI: A closed melonic graph Q has a unique melonic vertex pairing. 

P2: Consider cutting the two edges of some color {i} that are incident to the vertices of a melonic vertex pair. 
This splits the graph Q into two disjoint connected components (unless the edge of color {i} joined the 
vertices of the vertex pair directly). 

P3: Consider the k-bubbles of Q. These are also melonic and preserve the melonic vertex pairing defined by Q. 

P4: Consider a melonic vertex pair with vertices v and v, then v lies in a face of Q if and only if v lies in that 
face. In other words, they lie in the same D{D -\- 1)/2 faces. 

A. 1.3 1— dipole creation 

Let us remind ourselves that: 

Definition A. 3. A /c-dipole dk is a subset of Q comprising of two vertices v, v such that: 

— V and V share k edges of colors io, • • • , ^/c-i; 

— V and V lie in distinct {D -\- 1 — k)-huhh\es. In the arguments below, we say that the /c-dipole "separates" 

the vertices in the two bubbles. 

Lemma A. 4. Consider a melonic graph. Then, 1-dipole creation separates some melonic vertex pair. 

Proof. At the outset, let us note that closed {D + l)-colored graphs are 1-particle irreducible ioi D >1. Then, 
insert a 1-dipole of species {j} and suppose that it does not separate any melonic vertex pair. Thus, the two 
new vertices, v and v^ form a melonic vertex pair themselves. Therefore, by (P2), cutting the edges of color 
{z}, with i G {j}, emanating from v and v disconnects the graph into two connected components. But one 
of these cuts lies in each of the I^-bubbles of species {j} (separated by the 1-dipole). Thus, the I^-bubble is 
disconnected by a single cut, which contradicts the fact that it is 1 Pi. D 

Now, consider a melonic vertex pair vv. By property (P4), they both lie in the same D{D ^l)/2 faces, one of 
each species {j, k}. Now, let us examine these faces more closely, and in particular, their bounding edges: 

— £^J^ is the set of edges of color {j} that lie in the boundary of the face of species {j, k}. 

— SJ'^ — ri/eGl?! ^7^fe ^^ ^^^ ^^^ ^^ edges of color {j} that lie in the boundary of all such faces. 

Definition A. 5. Consider a generic (not necessarily melonic) graph that possesses a melonic vertex pairing. 
^0 ^dipole denotes all the edges that may be cut to create a 0-dipole, which separates the melonic vertex pair vv. 
^1 ^dipole i denotes all the edges that may be cut to create a 1-dipole of species {i}, which separates the melonic 
vertex pair vv. 

Lemma A. 6. For a melonic graph: 

cvv _ I I cvv cvv _ I I cvv fKQ\ 

^0 dipole ~ U i ' ^l-dipole,i ~ {J ^j * W^/ 

je{0,...,D} je(i} 
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Figure 17: A closed melonic graph, with melonic vertex pair vv and iPl melonic decorations. 

Proof. Consider a closed melonic graph in Figure 17, with melonic vertex pair vv. We have decomposed all 
melonic decorations into their iPl components, ^o^dipoie contains all, and only, those edges exterior to all iPl 
decorations. That ^o^dipoie cannot contain edges interior to some iPldecoration follows very simply from the 
irreducibility property (Q2). 

Similarly, ^ ^^dipoie,i contains those edges of color {j} (with j G {i}) exterior to all iPl decorations. D 

Finally, a result from [23] : 

Proposition A. 7. The set of rooted melonic graphs M is closed under 1-dipole creation and annihilation. 

A. 2 Identifying nlo core graphs 

At first glance, it is perhaps unsurprising that the graphs illustrated in Figure 7, being also the next simplest 
in terms of number of vertices, are the NLO core graphs. However, it emerges that it requires quite some care 
to prove this definitively. Closer inspection reveals that the graphs have degree: 

^^(^2 dipole) = ^-^^^{D - 2) . (60) 

Since the degree of ^2 dipole scales factorially with D, it is unclear at the outset that no other graph, perhaps 
less superficially obvious, sneaks in with a smaller degree to take the role of the NLO core graph. As it stands, 
the constraints on core graphs are not strong enough to rule out this possibility; one has only that they must 
be connected with p > 1. The aim of this subsection is to develop constraints that rule out this possibility. 

Lemma A. 8. Consider a {D -\- l)-colored core graph Q, at order p and D >A, with the following properties: 

- Q possesses exactly two melonic D-huhhles, say of species {zq} ci'^d {ii}. 

- Q possesses a planar jacket. 

Then, Q is the supermelon graph with p — 1 2-dipoles of species {iQ^h} inserted. 

Proof. Since ^ is a core graph, it contains a single I^-bubble for each color, labelled 'Tlr^ for z G {0, 1, . . . , D}. 
We shall denote the two melonic I^-bubbles by ^a )^ ^f?)' i"cspectively. Moreover, these I^-bubbles each 
contain all the vertices of Q. 

Since 7^/^ ^ and 7^/^ ^ are closed melonic graphs, by (PI), they each have unique melonic vertex pairing. Thus, 
a priori, a given vertex in Q has two melonic vertex pairings: one induced by 7^/^ ^ and another by 7^/^ y 
However, given the property (P3) of melonic vertex pairs, these vertex pairings coincide since both 7^.^ ^ and 

7^/^ N contain the {D — l)-bubbles of species (ioii)- Thus, even though Q is not melonic it has a melonic vertex 
pairing. 

Moreover, since Q possesses a planar jacket, we can draw it on the plane without crossings, such that the colored 
edges incident to each and every vertex are ordered according to some {D + l)-cycle r. 

We utilize the planar illustration of Q in Figure 18, where the cycle has the form r = (. . .jkl . . .). Consider a 
melonic vertex pair in Q and say that they do not share an edge of color k G {^o^^i}- Our argument has three 
subcases: 

i) Neither j nor / equal zq- ^^ '^u ), property (P2) ensures that cutting the edges of color k disconnects a 
subgraph from the rest of 7^/^ y Thus, there is a 1-dipole of species k in IZ/^ y Upon reinserting the lines 
of color zi, the planarity of the jacket ensures that this 1-dipole of species k persists into Q. 
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Figure 18: Killing 1-dipoles of species k. 

ii) Neither j nor / equal ii. The argument here is similar to the above with ii swapped for zq. 

^^0 {i'O — {^Oih}- In ^(i ), property (P2) ensures that cutting the edges of color k disconnects a subgraph 
from the rest of 1Zrj^ ^ As a result, the two distinguished edges of color k lie in the same face of color 
(ilk). Swapping the roles of io and zi, we can show that the two distinguished edges of color k lie in the 
same face of color {iok). Thus, back in ^, there is a 1-dipole of species k. 

Thus, all lines of color k ^ {^0,^1} directly connect melonic vertex pairs. In turn, this means that 7^/^ ^ is the 
I^-colored supermelon with p — 1 1-dipoles of species ^2 inserted (and vice versa). We illustrate 7^/^ ^ in Figure 
19. 




Figure 19: Distinguishing good from bad choices for the reinsertion of edges of color ii. 

Once again, the planarity of the jacket ensures that, upon reinserting the edges of color io, Q is the (Z) + 1)- 
colored supermelon with p — 1 2-dipoles of species {zq, ^1} inserted. 

D 

A corollary of this statement is the following: 

Corollary A. 9. Consider a {D -\- l)-colored core graph Q with D > A. If Q possesses three (or more) melonic 
D-bubbles and a planar jacket, then it is the supermelon. 



Now for the main result of this subsection, the identification of the NLO core graphs: 

Proposition 3.1. Consider the IID tensor model with D >3. The graphs S2 dipole o^^e the NLO core graphs. 

Proof. One uses an inductive argument on the number of colors. A /c-colored core graph is denoted by ^^"^Q. 
To start off, one shows that the statement holds true for D = 3. While Lemma A. 8 does not apply to this case, 
luckily it is simple to show directly. The following three statements hold in D = 3: i) Cc;('^^^^2 dipole) = 1, so it 
is certainly a NLO core graph (it is only sub-dominant by one power of N)] ii) ^^^^^2 dipole is the only core graph 
with p = 2; Hi) equation (19) implies that uj{^'^^Q) > p — 1, so given the second point, any other graph has a 
higher degree than ^^^^2 dipole- Thus, the NLO core graphs are as proposed: ^^^^nlo = ^^^^2 dipole- 



Say next that the statement holds true for the I^-colored model: *^^^Snlo = ^^■'G2-b 
i^(D-3). 



,. Thus w((^)eNLo) = 
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Now, say that despite this assumption, the statement does not hold for the {D + l)-colored model. Alas, the 
graphs ^^^"'^^52 DiPOLE, are trumped by some other graph ^^^"'^^ ^sneaky- In other words, their degrees satisfy the 
following inequality: 

U;((^ + ^)aSNEAKv) < C.((^ + l)a2-BIP0LE) = ^^^^ (^ " 2) • (61) 

A generic core graph of the {D + l)-colored model satisfies the equality: 

c.((^+i)a) = ^^y^(p - 1) + E-I^^'^^W) • (62) 

i 

where ^^^lZ(^i) are its I^-bubbles. In particular, *^^^^^Ssneaky obeys such a relation. They key now is to put a 
lower bound on both the number of vertices and the degrees of the I^-bubbles and force a contradiction. 

First of all, one knows also that p > 2, since the only core graph with p = 1 is the supermelon. Secondly, 
(^+^)^gj^EAKY possesses at least one planar jacket, else cc;(*^^+^^ ^sneaky) = ^j9j ^ D\/2. Moreover, appealing 
to Corollary A. 9, one can exclude all cases with three or more melonic D-bubbles. Thus, given that there are 
D + 1 D-bubbles in total, the configuration minimizing (62) has D — 1 oi them at NLO. But remember that 
i^-bubbles are i^-colored graphs, so the lower bound for the degree of ^^^^^Ssneaky is already: 

cu((^+^)asNEAKv) > ^^^^(p-l) + (D-l)^^^^(D-3)+c^r)%„))+cu((^)%,)) 
= ^^^{P + D-A)+ u;((^)7^(,,)) + c.((^)7^(,,)) . 

In order to satisfy the bound (61), one must set p = 2 and ijj{^^^lZ(^i^)) — uj{''^^1Z(^i^^) = 0. Thus, one falls into 
the case dealt with by Lemma A. 8 and a contradiction follows thereafter. D 

A. 3 Generating all NLO graphs 

From these core graphs, one can generate all graphs at NLO by performing arbitrary sequences of 1-dipole 
moves. At LO, such sequences resulted in the insertion of some set of elementary melons into the core graph. At 
NLO, while it is certainly true that most 1-dipole moves still result in the creation/annihilation of elementary 
melons, there is a subset that departs from this rule. The set of graphs S is certainly obtained by performing 
some sequence of 1-dipoles moves on the NLO core graph ^2 dipole- Our aim is to show now that these are all 
the graphs. 

With these properties at our disposal, we can proceed to analyse the NLO sector. 

Definition 3.2. One denotes by Si^n, the set of closed graphs derived from ^2 dipole,^ by inserting arbitrary 
combinations of n elementary melons. 

Proposition 3.3. The set of graphs: 

£>1 
n>0 

is closed under 1-dipole creation and annihilation. 

Proof. Every element of S^^n is built from the graph ^2 dipole,^ decorated by some combination of n melons. A 
generic graph in S^^n is drawn in Figure 20. 

1— dipole creation: From Definition A. 3, to create a 1-dipole of some species {z}, one must pick D edges of 
distinct colors taken from the set {i} that ensure the separation property. It will emerge that any 1-dipole 
involves D edges from a single partial subgraph of the type in Figure 21. 

To begin, the analysis of Section A. 2 shows that the graphs in Si^n have two melonic I^-bubbles (of species 
{io} and {^i}), which induce a unique melonic vertex pairing. ^^ One can easily identify this melonic vertex 
pairing in Figures 20 and 21. While the rooted melonic insertions have their melonic vertex pairs, those vertices 
exterior to any melonic decoration are in pairs of type AA, DD and so forth. 



■^^ For D = 3, there is little subtlety in that the NLO graphs have four melonic 3-bubbles. For all graphs in S^^n with i > 1, 
this poses no problem. However, examining the core graph ^2 dipole explicitly, one notices that it is symmetric and does not have 
a distinguished melonic vertex pairing. For this subset of graphs there are two different choices of pairing. This subtlety will be 
important later. 
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Figure 20: A generic element of 8^^^ where n counts 
the number of elementary melon insertions. 



Figure 21: A partial subgraph. 



In general, the creation of a 1-dipole in the graph corresponds to the creation of a 0- or 1-dipole its melonic 
/^-bubbles. By Lemma A. 4, a 1-dipole insertion separates some melonic vertex pair of a melonic I^-bubble 
(and a 0-dipole certainly does). Thus, a 1-dipole insertion separates some melonic vertex pair of the graph. 

Recall from Definition A. 5: ^i^dipoie i denotes the set of edges that may be cut to create a 1-dipole of species 
{%} in the graph, so as to separate the melonic vertex pair vv. 

One proceeds on a case by case basis. 

Case 1: Say one wishes to create 1-dipole of species {%} to separate a melonic vertex pair of type AA. 

The key is to examine the effect of the corresponding 0-/l-dipole on the melonic I^-bubbles. This restricts 
the elements of ^i^dipoie*- ^^^ finds that in all cases, ^i ^dipoie i contains those edges of species j (with 
j G {%}) in Figure 21 that are exterior to any melonic decoration. Any choice of D distinctly colored edges 
from this set sends the graph in S^^n to a graph in Si^n+i- 

Without loss of generality, let us assume that the melonic decoration between Bj and Bj is iPl, so that 
Bj and Bj form a melonic vertex pair. If it is not, then one decomposes this decoration into its iPl 
components and proceeds with one of the resulting vertex pairs. 



Case 2: Say one wishes to create a 1-dipole of species {%} to separate a melonic vertex pair of type BjB 



3' 



j G {^0, ^i}' Unless I^ = 3 and ^ = 1, ^i dipoie i contains only edges interior to the rooted melonic subgraph 
ABjBjA. Thus, we may apply Proposition A. 7 and one sees that 1-dipole creation sends the graph 
from Si^n to a graph in 5'^,n+i- The case D = 3 and ^ = 1 is special, as explained in footnote 14, since 
there are two possible pairings for the four vertices exterior to any melonic decoration. However, a 
pairing can always be chosen so that the assumption j G {zq, ^i} is valid and therefore we can avoid 
this case. 



J 



^o.^e{i^}: S^ 



pole i contains only edges interior to the rooted melonic subgraph ABqBqD. 
above, we may apply Proposition A. 7. 



As 



j = ii, i G {io}- As above, with the roles of io and ii swapped. 

j = io5 i = ii' As elements of ^i dip°oie i ' those edges of color {k} with k G {io} are interior to the 
melonic decoration between Bq and ^o- Those edges of color {i^} lie along the paths ABqBqD and 
DCqCqA. If the melonic decorations BqBq and CqCq are iPl, then the viable edges of color {io} are 
explicitly: ABq, BqD, DCq, CqA. If these melonic decorations are not iPl, then the edges of color 
{io} separating the iPl components must be added to the viable set. If one chooses to cut an edge 
of color {io} contained in the path ABoBoD^ then one sends the graph in Si^n to a graph in 5'^,n+i- 
However, if one chooses an edge of color {io} contained in the path DCoCoA^ then one sends the 
graph in S^^n to a graph in S^-^-i^n-i- This is the type of move illustrated in Fig 9. 

j = ii, i = io- As above, with the roles of io and ii swapped. 

Case 3: A 1-dipole to separate a melonic vertex pair nested at least once inside a melonic decoration. In 
this case, the viable edge set is contained entirely within the melonic insertion. Thus, we may apply 
Proposition A. 7 yet again. 

1— dipole annihilation: This involves picking an edge of color {i} such that after contraction, the number of 
I^-bubbles of species {i} is reduced by one. 
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Case 1: One picks a 1-dipole edge interior to some melonic insertion. Then, Proposition A. 7 holds and the 
graph in S^^n is sent to a graph in S^^n-i- 

Case 2: One picks a 1-dipole edge exterior to any melonic insertion. Then, the only 1-dipole edges are those 
or color {io} and {ii}. Consider the path A ^ Bq ^ Bq ^ D. There are two subcases: 

— The path from A to I) in Figure 21 is decorated by a melonic insertion. Then, the edge of color {zq} 

joining vertices A and Bq is a 1-dipole edge. Contracting this 1-dipole sends the graph in Si^n to a 
graph in Si^n-i- 

— The path from A to I) is undecorated. Then, the edge joining the vertices A and i^ is a 1-dipole edge, 

whose contraction sends the graph in S^^n to a graph in Si-i^n-\-i' This is the inverse of the move 
illustrated in Figure 9. 

D 

By cutting some single edge of the graphs in S, one generates the set of NLO connected 2-point graphs, denoted 
by T. These are defined and catalogued in Section 3.2. 

Proposition 3.4. For the IID model, the set of NLO connected 2-point graphs is: 

£>1 
m>0 
n>0 

Proof. This is a proof by inspection. One examines the graph in Si^n^ as drawn in Figure 20. For our purposes, 
there are two types edges in such a graph: 

— Cutting an edge interior to some melonic insertion sends a graph in S^^n to some graph in Ti^rn,n-mi where 

ra denotes the level of nesting (within the melonic insertion) of the cut edge; 

— Cutting an edge exterior to any melonic insertion sends a graph in 5*^^^ to some graph in T^,o,n- 
Moreover, taking a graph any graph in Ti^rn.n and joining its two open edges produces a graph in Si^n+m- D 
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